Reconstruction of Plane Circular Currents from Their Orthogonal Magnetic Field 
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The solution is presented to the problem what distribution of electric current in thin circular film 
provides a given distribution of normal (perpendicular to film) component of the current-induced 
magnetic field at film's surface. 



I. INTRODUCTION 

Plane electric current distributions and magnetic fields 
produced by them are of great interest in a number of ap- 
plications, for instance, when investigating magnetic flux 
pinning and critical states in thin superconductor films. 
At given current, one can calculate its field anywhere by 
direct integration. But in practice the inverse problem 
arises: how one could extract the current from measur- 
ing orthogonal (normal to plane) component of its field in 
some film's neighborhood, e.g. at its surface? Although 
this is standard problem of the classical theory of poten- 
tial, I did not find its solution in textbooks. Therefore 
I ventured to derive it independently, for special case of 
axial symmetry. 



II. THE PROBLEM 

Let a film with radius R and small thick ness D < < R 
(D/R -> 0) lies in the AF-plane (p = y/x 2 + y 2 < R , 
\z\ < D/2) and carries a circular current whose density 
integrated over the thickness is J(p) . We assume that 
in fact we know not J(p) itself but normal component 
of its magnetic field close to the film's surface, H z (p, z = 
±D/2 — * 0) . In case of superconductor film, this quan- 
tity represents density of vortices which pierce the film. 
From known H z (p,0) we want to obtain tangential ra- 
dial component of same field, H p (p,z = D/2 — > +0) = 
—H p (p,z = —D/2 — > —0) . If the latter is known, then 
the integral surface current density J(p) = J j<j>(p, z)dz 
(with iff, being azimuth current density) can be found 
with the help of relation 



4tt d f D/2 

— J(p) = ± 2H p (p, ±D 2) - — / H z (j>, 

C OP J -D/2 



, z)dz 



±2H p (p,±0) 



(1) 



which follows from the Maxwell equation curlff = 
Airj/c (written in the CGSE units). 

Outside the film, the field has potential character, 
H = —VU , with a potential U satisfying the Laplace 



equation AU = 0. With reference to circular geome- 
try, it is natural to consider the potential U using the 
spheroidal coordinates. Their definition and examples of 
there application can be found e.g. in . We use slightly 
modified version of spheroidal coordinates, in the form 



p 2 = R 2 {1 + u 2 ){l - v 2 ) , z 2 = R 2 u 2 v 2 
(0 < u 2 < oo , 0<v 2 < 1) , 



(2) 



at that the horizontal angle stays the third coordinate: 
x = pcos 4> , y — psin 4> . To write the Laplace operator, 
A , in these coordinates, it is convenient to introduce two 
operators 



Then A looks as 
1 



R z (u z + v z ) p 



i d 2 

~2 PlA/2 



(4) 



The film occupies the region p < R, z = , which cor- 
responds to u = , while the rest of the plane z = to 
v = . Of course, the potential must be anti-symmetric 
with respect to this plane, therefore, 



U(p > R, z = 0) = U(u, v = 0) = 



(5) 



At film's surface the potential has discontinuity, but its 
normal derivative is continuous. By the terms of our task, 



^-U(p<R,z) 



z=0 



H z0 (p) = H 



zQ 



Rv ou 



(Ry/l - « 2 ) 



u=0 



(0) 



where H z o(p) = H z (p < R,0) is a given function. For- 
mally, the problem AU = with boundary conditions 
and © is equivalent to the problem about electric po- 
tential of a charge distributed, with density H z0 5(z)/2tt , 
in the circular hole p < R cut out of ideally conducting 
plane z = 0. 

If the axial symmetry is assumed, then the Laplace 
equation reduces to 
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(A_ + A+)U = 



(7) 



2 



III. MAGNETIC POTENTIAL 

The operator A + in is nothing but generating op- 
erator for the complete set of orthogonal (at interval 
—1<v<1) Legendre polynomials P n (v) (see e.g. @): 



A+(v)P n {v) = -X„P n (v) , A„ = n(n + 1) 



(8) 



Therefore, particular solutions to Eq[7| have the form 
Pn(v)Q n (u) , where functions Q n {u) satisfy the equa- 
tion A_(u)<2„(u) = X n Q n (u) . It is easy to see that 
A_(w) = — A + (i>) . Consequently, Q n (u) oc P n (iu) , 
with P n (u) also being solutions of Eq|SJ A + (u)P n (u) = 
—X n P n (u). That must be chosen be the second kind 
(non-polynomial) eigenfunctions which have zero asymp- 
totic at infinity, because U(u,v) should tend to zero at 
z — > oo. Hence, according to general theory of classical 
special functions Q, Qn(u) can be represented by 



!2k — 



1 uP 2 k(s)ds 



>2k+l = 



sP 



2fe+l 



{s)ds 



s 2 + u 2 



(9) 



In view of the boundary condition (JSJ , as well as (0 , 
the odd polynomials with n = 2k + 1 (k = 0, 1, ..) only 
contribute to complete solution of Eq[7| By this reason, 
we can write 



with coefficients C k to be determined from condition 10 . 
Notice 2] that P 2k+ i(v) form complete orthogonal set 
in the half-interval < v < 1 . Under standard classical 
definition of Lcgcndrc polynomials, 

P 2k+1 (s)P 2m+1 (s)ds = 5 km /(4k + 3) , (11) 

that is 4fc + 3 in (|10fl serves as normalizing multiplier. 

The result of subsequent proper calculations (detailed 
in Appendix in Sec.l) reads 

Ck = (-l) k S k £ P 2k+1 (yi--^j H zQ (r) r -^ (12) 

Here the quantities 

S k = £(1/2, k + l)/2 = 2 k kl/(2k + 1)!! (13) 

are introduced. At upper side of film's surface, formulas 
(fTUf) and (O, as combined with © and ©, yield (see 
Appendix 2) 



U(p,z = +0) = G R (p,r)H z0 (r)rdr , (14) 



U(u, v)=J2 ( 4fc + 3 ) Cfc -ftfc+i (») Q2fc+i (u) , (10) where the G reen function G R (p,r) is presented by 



k=0 



R 2 



R 2 



(15) 



7r max(p, r) 



A' 



min(p, r) 
max(p, r) 



— F ( arcsin 



max(p, r) min(p, r) 
i? ' max(p, r) 



(16) 



7r (/? + r) 



AT 



2^p¥ 



2^ - arcsin — arcsin — , — 

p+r j \2 R 2 R' p + r 



(17) 



where F(</>, A:) is the first-kind elliptic integral and 
K(k) — F{n/2 : k) the first-kind complete elliptic inte- 
gral. The way of summation of the series l|15ll is described 
in Sec. 3 of Appendix. 



IV. CURRENT 

Provided the normal field at surface of (infinitely) thin 
film is known, the current can be restored by means of 



and ltl4T> : 



2tt , N 
-J(P) 



H p o(p) = --i L l G R {p,r) H z0 (r)rdr , 
a P Jo 

(18) 

-0) . Integration of fT£|) 



where H p0 (p) = H p (p,z = 
yields the mean current: 



2tt 1 
~~c R 



Jd P =- 



2 . r . 
— arcsin — I H z q (r) dr 



R 



(19) 



It is easy to calculate also the magnetic moment, M , of 
the current 0: 



M = | J Q H r H z0 (r) dr (20) 

Thus both these characteristics eliminate contributions 
from the film's edge, i.e. H z0 (p ss R) . 

The expansion when substituted to l|18[) . deter- 

mines "natural modes" of the surface field: 



H z0 = h ±k (p) 



P 2k+1 (^l-(p/R) 2 ) 

Vi - (p/^) 2 



pO 



h\\k(p) = ~S k 



2 c2 pP^Wl-ip/Rf 



(21) 



(22) 



Clearly, h± k (p) is 2fc-order polynomial, and h± k (0) = 1 
(because P n (l) = 1), while h^ k (p) is polynomial of 
the same order multiplied by factor pj \J R 2 — p 2 and 
possesses square-root divergency at film's edge. Both 
h±k(p) and h\\ k (p) have k zeros at < p < R . In 
particular, h\\ (p) — 2p/n\J R 2 — p 2 determines the cur- 
rent distribution (see e.g. 4i|) which creates constant 
unit-value field at film's surface. Of course, one can form 
linear combinations of these modes without current's di- 
vergency at the edge. For example, the two lowest modes 
( k = 0, 1 ) give 

3 



H z0 = 1 



^ H, 



(23) 



with x = p/R . 

The approach mentioned in Appendix 3 allows to di- 
vide (|18f) into singular and regular parts: 



H P0 (P) = HT 3 (P) + HT(P) , 



h s p 7 9 (p) 



2p 



p* Jo 



R H z0 (r)rdr 
R.^R 2 - r 2 



p r\r dH z0 (r) 



dr 



dr 



In the latter formula, 



(24) 



(25) 



(26) 



(27) 



_ 2 [K(k) - F(<j>, k) - E{k) + Efa k)} 

1 \ x i V) — ■ i \ 

it mm(x, y) 

with E(k) and E(<fi, k) being the second-kind elliptic 
integrals whose module of ellipticity k and phase <f> are 
expressed by 

k = min(a:, y)/ max(cc, y) , <fi = arcsinmax(x, y) (28) 

Clearly, T(x, y) has logarithmic peak at x = y (see 
Fig.l). What is important, T(l, y) = T(x, 1) = , there- 
fore H r p Q 9 (p — > R) is zero (or at least finite value). In 

opposite, Hp™ 9 (p — ► R) diverges except the cases when 
the integral in (|25|) turns into zero (e.g. in the example 
(J2SJ). In such case, (I18fl becomes integral relation l|26() 
between current and radial derivative of normal field (in 
place of differential relation Anj/c = —dH z /dp which 
would work under cylindrical geometry). 




FIG. 1: The kernel T(x,y)y of the integral relating the 
current to the field gradient, as function of y at x = 
0.1, 0.3, 0.5, 0.7 and 0.9. 



V. CONCLUSION 

To resume, we found how to determine an electric cur- 
rent distribution in thin circular film from normal compo- 
nent of the current-induced field measured in the vicinity 
of film's surface (very similar task about infinite strip film 
must be considered separately). The results can be ap- 
plied, in particular, to numeric modeling of magnetic flux 
penetration into superconducting film. 

I am grateful to Dr. Yu. Medvedev and 
Dr. V. Khokhlov for useful comments. 



APPENDIX 



1. According to ©, HHJ) and (JTTJ, 

Cfe = -IQWit+or 1 x 

x R £ s P 2k +i (s) H z0 (RVT^) ds , 



(A.l) 



where the upper half-space z > is considered, and the 
prime means differentiation: 

Q'ik+i (+°) = lim u^+o dQ 2k+l {u)/du (A.2) 
From @ it follows that at 

» , , /, m, I 1 u2p 2k+l(s)ds 
hk+l(U) = Q2k+1{0) - 



_l s(u 2 + s 2 ) 

Here in the limit u — > +0 factor u/(u 2 + s 2 ) works as 
ttS(s) , therefore 

Q2 fe+ i(+0) = -7T lim P 2k+1 (s)/s = -nPi k+1 (0) (A.3) 

The right-hand sides of these equalities can be easily 
found from the generating function of Legendre polyno- 
mials 0, 



t n P n (v) = l/y/l + t 2 - 2tv 



(A.4) 



4 



Its differentiation by v at point v = and Taylor ex- 
pansion over t show that 



^2 fc+ i(0) = (-l) k /S k , 



(A.5) 



with Sk defined by JT^Jl, and P^fc+ i ( s ) = dP2k+i{s) / ds . 
Combining (| A. 1|) . (|A.3(1 and (|A.5(1 and performing the 
obvious change of variables in the integral in (|A.1(I . we 
obtain (fT2")) . 

2. According to (|1UI) and (J2J, to express the po- 
tential at upper film's side, U(p,z — > +0) = C/(w — > 
+0, i> = yl — p 2 / R 2 ) , we should calculate the quanti- 
ties Q2fe+i(0) determined by This can be made ei- 
ther with the help of the recurrent relations between Leg- 
endre polynomials or with the help of i|A.4|) . Direct 
integration of (|A.4() . after dividing by v and replacing t 
by it , yields 



- T m 2k+l I 1 P2k+l{v) dv = aresin I 



fe=0 



(A.6) 



From another hand, let us pay attention to that the co- 
efficients <|13[1 can be represented as 

S k = f* /2 sm 2k+1 6d8 , (A.7) 
therefore their generating function coincides with i|A.6|) : 



fc=0 



i - t 2 sin 2 e» VT^T 2 



Comparing (|A.8|I and 1A.6|) . and taking into account @ 
and l|13[) . we conclude that 



2fc +i(o) 



1 P 2 fc+i(s) 



'-1 s 

3. Let us consider the series 



ds = 2(-l) k S k (A.9) 



is based on the remarkable formula [!|: 



I Re 1/vf 



(A.ll) 



y2 _ ^2 _|_ 2a;y2; 



where —1 < x, y,z < 1 and the symbol Re means tak- 
ing real part, that is the right-hand side is zero when the 
subradical expression is negative (in |3] this formula is 
presented in slightly different notations). 

Now, divide both sides of (|A.11(I by factor 2z and 
integrate over z from —1 to 1 (in the sense of principal 
value). On the right-hand side we have exactly calculable 
integral. On the left, apply the relation (| A. 9|> (clearly, the 
odd terms only survive, with indices n = 2k + 1). The 
result is 



= Re sign(a;y) / y 7 ^ 2 + y 2 — 1 

(A.12) 

Next, multiply together g(x, s) and g(y, s) and in- 
tegrate the product over < s < 1 . On the left we 
should apply the orthogonality and normalization rela- 
tions Q11JI. On the right, we arrive to a standard elliptic 
integral, eventually obtaining the equality 



a(x,y) 



sign(x;/) [K(k) - F((f>, k)] 
y 7 ! — mm(x 2 ,y 2 ) 



(A.13) 



where K(k) and F((f>,k) are standardly designated first 
kind elliptic integrals (complete and incomplete, respec 
tively) whose module of ellipticity k and phase 
expressed by 



are 



arcsin y 7 ! — min(a; 2 , y 2 ) 



k — y[l — max(cc 2 , y 2 )] / [1 — min(x 2 , y 2 )] 



oo 

a{x, y) =J2 (4fc + 3) S 2 k P 2k+ i(x)P 2k+1 (y) , (A. 10) and function ^y) is de fined b y (£3®. C hanging here 



fc=0 



which turns into H15fl after the evident change of variables 
(and adding constant multiplier). One of ways to sum it 



x and y by y/l — p 2 /R 2 and — r 2 / R 2 , respectively, 
and multiplying the result by 2/ nR , we come to the 
relations lfT3li - lfTCli . The conversion of (|To)l into ltT7|l is 
based on properties of elliptic integrals (see e.g. pjlfl)- 
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